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Abstract 

This paper is on beamforming in wireless relay networks with perfect channel information 
at relays, the receiver, and the transmitter if there is a direct link between the transmitter and 

C^^ , receiver. It is assumed that every node in the network has its own power constraint. A two-step 

amplify-and-forward protocol is used, in which the transmitter and relays not only use match 

Tj- ' filters to form a beam at the receiver but also adaptively adjust their transmit powers according 

o. 

00 , to the channel strength information. For a network with any number of relays and no direct link, 

the optimal power control is solved analytically. The complexity of finding the exact solution 

/\ ■ is linear in the number of relays. Our results show that the transmitter should always use its 

maximal power and the optimal power used at a relay is not a binary function. It can take any 
value between zero and its maximum transmit power. Also, surprisingly, this value depends on 
the quality of all other channels in addition to the relay's own channels. Despite this coupling 
fact, distributive strategies ai^e proposed in which, with the aid of a low-rate broadcast from the 
receiver, a relay needs only its own channel information to implement the optimal power control. 
Simulated performance shows that network beamforming achieves the maximal diversity and 
outperforms other existing schemes. 
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Then, beamforming in networks with a direct link are considered. We show that when the 
direct link exists during the first step only, the optimal power control at the transmitter and 
relays is the same as that of networks with no direct link. For networks with a direct link during 
the second step only and both steps, recursive numerical algorithms are proposed to solve the 
power control problem. Simulation shows that by adjusting the transmitter and relays' powers 
adaptively, network performance is significantly improved. 

1 Introduction 

It is well-known that due to the fading effect, the transmission over wireless channels suffers from 
severe attenuation in signal strength. Performance of wireless communication is much worse than 
that of wired communication. For the simplest point-to-point communication system, which is com- 
posed of one transmitter and one receiver only, the use of multiple antennas can improve the capacity 
and reliability. Space-time coding and beamforming are among the most successful techniques de- 
veloped for multiple- antenna systems during the last decades [1,2]. However, in many situations, 
due to the limited size and processing power, it is not practical for some users, especially small 
wireless mobile devices, to implement multiple antennas. Thus, recently, wireless network commu- 
nication is attracting more and more attention. A large amount of effort has been given to improve 
the communication by having different users in a network cooperate. This improvement is conven- 
tionally addressed as cooperative diversity and the techniques cooperative schemes. 

Many cooperative schemes have been proposed in literature [3-21]. Some assume channel in- 
formation at the receiver but not the transmitter and relays, for example, the noncoherent amplify- 
and-forward protocol in [8,9] and distributed space-time coding in [10]. Some assume channel 
information at the receiving side of each transmission, for example, the decode-and-forward proto- 
col in [8, 12] and the coded-cooperation in [13]. Some assume no channel information at any node, 
for example, the differential transmission schemes proposed independently in [14-16]. The coher- 
ent amplify-and-forward scheme in [9, 1 1] assumes full channel information at both relays and the 
receiver. But only channel direction information is used at relays. In all these cooperative schemes. 



the relays always cooperate on their highest powers. None of the above pioneer work allow relays 
to adjust their transmit powers adaptively according to channel magnitude information, and this is 
exactly the concern of this paper. 

There have been several papers on relay networks with adaptive power control. In [22, 23], 
outage capacity of networks with a single relay and perfect channel information at all nodes were 
analyzed. Both work assume a total power constraint on the relay and the transmitter. A decode-and- 
forward protocol is used at the relay, which results in a binary power allocation between the relay 
and the transmitter. In [24], performance of networks with multiple amplify-and-forward relays and 
an aggregate power constraint was analyzed. A distributive scheme for the optimal power allocation 
is proposed, in which each relay only needs to know its own channels and a real number that can 
be broadcasted by the receiver. Another related work on networks with one and two amplify-and- 
forward relays can be found in [25]. In [26], outage minimization of single -relay networks with 
limited channel-information feedback is performed. It is assumed that there is a long-term power 
constraint on the total power of the transmitter and the relay. In this paper, we consider networks 
with a general number of amplify-and-forward relays and we assume a separate power constraint on 
each relay and the transmitter. Due to the difference in the power assumptions, compared to [24], 
analysis of this new model is more difficult and totally different results are obtained. 

For multiple-antenna systems, when there is no channel information at the transmitter, space- 
time coding can achieve full diversity [1]. If the transmitter has perfect or partial channel infor- 
mation, performance can be further improved through beamforming since it takes advantage of 
the channel information (both direction and strength) at the transmit side to obtain higher receive 
SNR [2] . With perfect channel information or high quality channel information feedback from the 
receiver at the transmitter, one-dimensional beamforming is proved optimal [2,27,28]. The more 
practical multiple- antenna systems with partial channel information at the transmitter, channel statis- 
tics or quantized instantaneous channel information, are also analyzed extensively [29-33]. In many 
situations, appropriate combination of beamforming and space-time coding outperforms either one 
of the two schemes alone [34-37]. In this paper, we will see similar performance improvement in 
networks using network beamforming over distributed space-time coding and other existing schemes 



such as best-relay selection and coherent amplify-and-forward. 

We consider networks with one pair of transmitter and receiver but multiple relays. The receiver 
knows all channels and every relay knows its own channels perfectly. In networks with a direct 
link (DL) between the transmitter and the receiver, we also assume that the transmitter knows the 
DL fully. A two-step amplify-and-forward protocol is used, where in the first step, the transmitter 
sends information and in the second step, the transmitter and relays, if there is a DL, transmit. We 
first solve the power control problem for networks with no DL analytically. The exact solution can 
be obtained with a complexity that is linear in the number of relays. Then, to perform network 
beamforming, we propose two distributive strategies in which a relay needs only its own channel 
information and a low-rate broadcast from the receiver. Simulation shows that the optimal power 
control or network beamforming outperforms other existing schemes. We then consider networks 
with a DL during the first transmission step, the second transmission step, and both. For the first 
case, the power control problem is proved to be the same as the one in networks without the DL. For 
the other two cases, recursive numerical algorithms are provided. Simulation shows that they have 
much better performance compared to networks without power control. We should clarify that only 
amplify-and-forward is considered here. For decode-and-forward, the result may be different and it 
depends on the details of the coding schemes. 

The paper is organized as follows. In the next section, the relay network model and the main 
problem are introduced. Section[3]works on the power control problem in relay networks with no DL 
and Section |4] considers networks with a DL. Section [5] contains the conclusion and several future 
directions. 

2 Wireless Relay Network Model and Problem Statement 

Consider a relay network with one transmit-and-receive pair and R relays as depicted in Fig. \T\ 
Every relay has only one single antenna which can be used for both transmission and reception. 
Denote the channel from the transmitter to the zth relay as /j and the channel from the zth relay 
to the receiver as gi. If the DL between the transmitter and the receiver exists, we denote it as 



/o- We assume that the transmitter knows /o, the ith relay knows its own channels fi and gi, and 
the receiver knows all channels /o, /i, . . . , /r and gi, . . . , g^. The channels can have both fading 
and path-loss effects. Actually, our results are valid for any channel statistics. We assume that for 
each transmission, the powers used at the transmitter and the zth relay are no larger than Pq and Pj, 
respectively. Note that in this paper, only short-term power constraint is considered, that is, there is 
an upper bound on the average transmit power of each node for each transmission. A node cannot 
save its power to favor transmissions with better channel realizations. 

We use a two-step amplify-and-forward protocol. During the first step, the transmitter sends 
ttoV^s. The information symbol s is selected randomly from the codebook S. If we normalize it 
as E|sp = 1, the average power used at the transmitter is o^Pq. The ith relay and the receiver, if a 
DL exists during the first step, receive 

Ti = ao^/PofiS + Vi and xi = aov^/o-5 + w^i, (1) 

respectively. Vi and W\ are the noises at the ith relay and the receiver at Step 1 . We assume that they 
are CA/'(0, 1). During the second step, the transmitter sends /SoV^e-^^^s, if a DL exists during this 
step. At the same time, the ith relay sends 



The average transmit power of the ith relay can be calculated to be afPi. If we assume that /o keeps 
constant for the two steps, the receiver gets 
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W2 is the noise at the receiver at Step 2, which is also assumed to be CJ\f{0, 1). Note that if the 
transmitter sends during both steps, we assume that the total average power it uses is no larger than 
Pq. With this, the total average power in transmitting one symbol is no larger than Xli=o -^«- Clearly, 
the coefficients ao, ai, ... ,aji are introduced in the model for power control. The power constraints 
at the transmitter and relays require that Oq + /9q < 1 and < Oj < 1. 



Our network beamforming design is thus the design of 9o, 6i, ■ ■ ■ ,6r and ao, /5o, ai, ■ ■ ■ , an, 
such that the error rate of the network is the smallest. This is equivalent to maximize the receive 
SNR, or the total receive SNR of both branches if a DL exists during the first step. From (O, we 
can easily prove that an optimal choice of the angles are Oq = — arg /o and 9i = — (arg /j + arg^fj). 
That is, match filters should be used at relays and the transmitter during the second step to cancel 
the phases of their channels and form a beam at the receiver. We thus have 






What is left is the optimal power control, i.e., the choice of ao, /^o, «i, • • • , «r- This is also the main 
contribution of our work. 

3 Optimal Relay Power Control 

In this section, we investigate the optimal adaptive power control at the transmitter and relays in net- 
works without a DL. Section [STI presents the analytical power control result. Section [3^ comments 
on the result and gives distributive schemes for the optimal power control. Section 13.31 provides 
simulated performance. 

3.1 Analytical Result 

With no DL, we have (3o = and xi = 0. From ([3]), the receive SNR can be calculated to be 

^, JD ( \^R ai\fi9i\VP~i \ 
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It is an increasing function of uq. Therefore, the transmitter should always use its maximal power, 
i.e., ftQ = 1. The receive SNR is thus: 
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Before going into details of the SNR optimization, we first introduce some notation to help 
the presentation. (■, ■) indicates the inner product. || ■ || indicates the 2-norm. P indicates the 
probability. Oj denotes the ith coordinate of vector a and aj^ ... j^, denotes the fc-dimensional vector 
where -^ represents the transpose. If a, b are two /^-dimensional vectors, 
1, . . . , R. Or is the i?-dimensional vector with all zero entries. 
Denote the set 0/? ^ y ^ a or equivalently, < yi < ai for i = 1, . . . , R, as A. For 1 < k < R— 1, 
denote the set 0^ ^ yiu...,ik ^ ^i,-,ik as An,...,**,, where {zi, . . . , 4} is a /c-subset of {1, ... , R}. 
Define 



Gj-j^ ■ ■ ■ flj^ 

a ^ b means a, < 6, for all i 
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and A = diag{a}, 



where diag{a} indicates the diagonal matrix whose ith diagonal entry is cii. With the transformation 
y = Ax, or equivalently, x = A^^y, we have 



SNR = Po: 



(b,x)^ 



lAxI 



Pn 



(c,y)' 



where 



A-^h 



\/l+l/iPJ"o 
l9i|v^ 
















l/il 



The receive SNR optimization problem is thus equivalent to 

(c,y)^ 



max- 



s.t. y G A. 



(4) 



y l + ||y||2 

The difficulty of the problem lies in the shape of the feasible set. If y is constrained on a hypersphere 
that is, ||y|| = r, the solution is obvious at least geometrically. Given that ||y|| = r. 
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where Lp is the angle between c and y. The optimal solution should be the vector which has the 
smallest angle with c. Thus, we decompose dH) as 
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max- 

r 



I max (c,y) | s.t. y G A and < r < ||a||. (5) 

Since P(aj > 0) = 1 and P(cj > 0) = 1, we assume that aj > and q > 0. Define 

(6) 






for i = 1, . . . , R and, for the sake of presentation, define (f)R+i = 0. Order (pj as 



Vi > 0r2 > ■■• > (Prn > <Prn+i- 



(ti, r2, . . . , tr, T/j+i) is thus an ordering of (1, 2, . . . , i?, i? + 1) and r^j+i = i? + 1. Define 



ro=0, 



ri= 



^2 = 
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Since 0r _i > 0t , we have rj_i < r^ for j = !,...,/?. Thus, the feasible interval of the radius, 
[0, ||a||], can be decomposed into the following R intervals: 

[0> l|a||] = [ro, ri] U [ri, rg] U ■ ■ ■ U [r_R_2, r_R_i] U [r^j-i, r^]. 
We denote Fj = [rj, rj+i] for z = 0, . . . , i? — 1. Thus, ([5]) is equivalent to 

max max ( max (c, y) 

;=i,...,ij r-er, 1 + r^ V llyli='-er,,yeA 



We have decomposed the optimization problem into R subproblems. We now work on the ith 
subproblem: 



max max 

reFi 1 + r^ \ ||y||=r-Gri,yeA 



(c,y) 



Denote the solution of the inner optimization problem, 



max (c,y), 
llylNreri.yGA 



(8) 



(9) 



as z*^*\ We have the following two lemmas. 



Lemma 1. Zj = ajforj = ri, . . . , Tj. 
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Proof. We prove this lemma by contradiction. Assume that Zj < aj for some j G {ri, . . . , r^}. We 
first show that there exists an / G {rj+i, . . . , tr} such that -^— < -^. Assume that -^— > s^^ for all 
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Cj Ci 



m G {ri+i, . . . , tr}. We have zi^' < Cm^ < Cmj- = Cm4>j ^ Thus, 
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because of d?]) 



This contradicts ||z^*^|| G Fj. thus, there exists an / G {tj+i, . . . , tr} such that -^ < -^. 



Define another vector z' as z' = 2;, + 6, z\ 
m ^ i,l, where 

< 5 < min <! 2c, ( 1 + 41 V^^ - ^ 1 , 



.(0 



26 zf'^ - 52, and z' = z^l^ for 
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(i) 



.(») 



r(') ,(0 



Since we have assumed that z, < aj and have just proved that -jr- < -j—, such 6 is achievable. To 
contradict the assumption that z*^*) is the optimal, it is enough to prove the following two items: 



1. z' is a feasible point: ||z'|| = r and z' E A, 

2. (c,z«)<(c,z'). 

From the definition of z', we have 






)||2 ^^2_ 



Since < 6 < aj - zf\ we have < z^- < aj. Also, since < 6 < J (zf^) + [4') - zf , we 



can easily prove ,ha. - = ^(.«)' - 2fef - «^ > and .; < .,"> < a,. Thn., z' e A. The fir., 
item has been proved. For the second item, since 5 < 2cj \\ + -^ \ \—^ ^ J , we have 



Cizf + Cjzf - [Ciz\ + CjZ\) < 

(c,z«)<(c,z'). 



D 



Lemma 2. zf = ^^^^^^c.forj = r,+i, . . . , tr. 

Proof. From LemmafU zj = aj for j = ri, . . . , Tj. Thus, ^ can be written as 

i 

y =r-eri,yGA ^ — ' 
m=l 

i 
= Yl ^-- + "^ aX (Cr,+i,...,r«, yr,+i,...,rH)- 

m=l llyTi + l.-'-^flll^V '' "^m^l^Tm" 

''<=ri,yTi_,_i,...,Tjj6AT-^_l_-^,..._T-^ 



/ 2 v~^i 2 

Define A = ^L" "'";7" - It is obvious that (c^,+,,...,^^, y^,^,,...,^^) < (c^^^,,...,^^, Ac^,^,„„,^^) for all 



\'yTi+i,...,TR\\ = yr"^ — Z]m=i ^Tm- ^ other words, to maximize the inner product, yr,+i,...,Tfl should 
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have the same direction as c^^^^^ ,,, -r^. Thus, we only need to show that this direction is feasible for 
r G Fj. This is equivalent to show that Ac^-^^ ... ^-^ G A^-^^ ... ^-^ for any r G Fj. We can easily prove 
that 

r G F,- <^ A G Qi, 



where f2, 



K^K. 



for i = 0, . . . , i? — 1. Thus, for any r G Fj and j = Tj+i, . . . , t/?, we have 



< \cj < (p^^l^Cj < (pj ^Cj = ttj. Hence, Ac,,^,,...,,^ G Ar,+i,...,r«- □ 

Combining Lemma [Hand Lemma [2l we have 

and thus 

i 

max (c,y) = ^6,„ + A||c,,^,,...,,^f. (11) 

m=l 

We have solved the inner optimization of Subproblem i. The solution of the R subproblems can 
thus be obtained. 

Lemma 3. For i = 1, . . . ,R, define 

The solution of Subproblem is y^^^^ = (fi'^^c. The solution of Subproblem ifor i = l,...,R — lis 
y*^*-* that is defined as 

{ttj j = Ti,...,Ti 

mm |Ai, 0^^J_^ I Cj J = Ti+i, . . . , tr. 

Proof From (fTT)) . Subproblem i is equivalent to the following 1 -dimensional optimization problem: 

max — ^ ■■ — . (13) 

When i = 0, (1131) is equivalent to max^eQ^ il\cp\2 - Since j^ln^p^a is an increasing function of 
A, its maximum is at A 



n 
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For i = 1, . . . , R— 1, Define 

2 



/^V* h -J- \\n l|2X 

^^^ '^ 1 , V^^ 9 ill II 9 \ 

We have, 

d^i 2 (^Em=l ^r,„ + \\Cr,+,,...,rnf\) || C^,+i,...,r« || 



i + E<-E^-^ 



12^2' ^ m=l m=l 



Thus, ^ > if A < Ai and || < if A > A^. So, if Aj < 0~/^, the optimal solution is reached at 
A = Aj. Otherwise, the optimal solution is reached at A = (p~^ . From (flOl) . Subproblem i is solved 
at y('^ as defined in ^. D 

Now, we can work on the relay power control problem presented in dH). 

Theorem 1. Define x^*) as 

{1 j = ri,...,ri 

(14) 
Ai0j j = Ti+i,...,TR 

The solution of the SNR optimization is x'^*"\ where iq is the smallest i such that \i < (j^r^^- 



Proof. First, since 0r+i = 0, we have A/j < tp"^ ^ = 0^+i = oo. Thus, ?o exists. Also, since 
Ajo < 0~^ ^, and 0^-^ decreases with j, we have x^-°' < 1 for j = Tjo+i, . . . ,tr. This means that 
x(*o) is in the feasible region of the optimization problem. 
Denote 

Note that ||y*^°) || = ri. Since ri G Fi, y*^") is also a feasible point of Subproblem 1. Thus, 7] (y*^°^) < 
r] (y*-^') due the optimality of y*^^^ in Subproblem 1. This means that there is no need to consider 
Subproblem 0. For z = 1, . . . , _R — 2, if Aj > 0"^ , 
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and 
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N 



-2 \\n II _L \ ' 1/7 |2 



'-^T,- I 1 Tro ~r / f^T„- ' ' 



^r,+ i H^n+i,...,rfl|| ^/^l"r,l - ' i+1 



E 



i=i 
Since rj+i G Fj+i, yW is a feasible point of Subproblem i + 1. Thus, 7] (y'-*-') < ?7 (y*^*"'"^'') due the 

optimality of y(*+^) in Subproblem i + I. This means that there is no need to consider Subproblem 

z + 1. Thus, we only need to check those y(*)'s with A^ < 0^^ , and find the one that results in 

the largest receive SNR. From the definition in (fT4l) . this is the same as to check those x^'^'s with 

Now, we prove that Aj+i < ^"^^ if -^j < ^n^i- First, from Aj < 0^-^^, we have 

Ei 1 ^n+i' 

m=l "Tm 

Since 7-^^^ = 0_/^^ , we can prove easily that 

Thus, we only need to check those x*^*) 's for io < i < R and find the one causing the largest receive 



SNR. From previous discussion, io > I- 

Define SNRi = -^^. 
proof of Lemma 3, we have 



Define SNRi = j^^^- Now, we prove that SNRi > SNRi+i for Iq < i < R. From the 



2 



(\^^ h -I- llr- l|2\ ^ 

SNRi = ^ : ^— 

1 I Y^« ^2 _|_ llf, ||2\2 



V- ^2 (^™=1^ 



77T.=2+1 



1 + Em=l «L 



--SNRi+i ^ -^ ^ , , ^i r ~ 1 , V^»+l 



a. 



,2 



1 + E™=i «L 

Thus, the optimal power control vector that maximizes the receive SNR is x^*''^ D 
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3.2 Discussion 

It is natural to expect the power control at relays to undergo an on-or-off scenario: a relay uses its 
maximum power if its channels are good enough and otherwise not to cooperate at all. Our result 
shows otherwise. The optimal power used at a relay can be any value between and its maximal 
power. In many situations, a relay should use partial of its power, whose value is determined not 
only by its own channels but all others' as well. This is because every relay has two effects on the 
transmission. For one, it helps the transmission by forwarding the information, while for the other, 
it harms the transmission by forwarding noise as well. Its transmit power has a non-linear effect on 
the powers of both the signal and the noise, which makes the optimization solution not an on-or-off 
one, not a decoupled one, and, in general, not even a differentiable function of channel coefficients. 
As shown in Theorem [Hand Lemma [3l the fraction of power used at relay j satisfies Oj = 1 for 
j = Ti, . . . , Tjo and aj = Xifj(l)j for j = Tjq+i, . . . ,tr. Thus, the zq relays whose 0's are the largest 
use their maximal powers. Since iQ>l, there is at least one relay that uses its maximum power. This 
tells us that the relay with the largest always uses its maximal power. The remaining R — io relays 
whose 0's are smaller only use parts of their powers. For j = Tjq+i , . . . ,tji, the power used at the jth 



relay is a^^P, = Xlcf'Pj = A^J/,/^?,^ (1 + |/,|2Po), which is proportional to |/,/(?,f (1 + |/,f P, 



o; 



since Aj^, is a constant for each channel realization. Although Pj does not appear explicitly in the 
formula, it affects the decision of whether the jth relay should use its maximal power. Actually, in 
determining whether a relay should use its maximal power, not only do the channel coefficients and 
power constraint at this relay account, but also all other channel coefficients and power constraints. 
The power constraint of the transmitter, Pq, plays a roll as well. 

Due to these special properties of the optimal power control solution, it can be implemented 
distributively with each relay knowing only its own channel information. In the following, we 
propose two distributed strategies. One is for networks with a small number of relays, and the other 
is more economical in networks with a large number of relays. 

The receiver, which knows all channels, can solve the power control problem. When the number 
of relays, R, is small, the receiver broadcasts the indexes of the relays that use their full powers and 
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the coefficient Xi^. If relay j liears its own index from the receiver, it will use its maximal power 
to transmit during the second step. Otherwise, it will use power Xjjfj/gjl"^ (1 + |/,pPo)- The bits 
needed for the feedback is 

io\ogR + Bi < RlogR + Bi, 

where zq is the number of relays that use their maximal powers and Bi is the number of bits needed 
in broadcasting the real number Ajg. Instead, the receiver can also broadcast two real numbers: Ajg 
and a real number d that satisfies (pn > d > (p^ +i- Relay j calculates its own (pj. If cpj > d, 
relay j uses its maximal power. Otherwise, it uses power XfJ/j/gj]"^ (1 + \fj\'^Po). The number of 
bits needed for the feedback is 2Bi. Thus, when R is large, this strategy needs less bits of feedback 
compared to the first one. 

Networks with an aggregate power constraint P on relays were analyzed in [24]. In this case, 
with the same notation in Section [STl Pj = P and Xli=i '^i — l- The optimal solution is 



«i 






R \fmgm\Hl + \fm\^P0) 



E 



^"=1 (l/m|2Po + |9™|2P+l) 



aj is a function of its own channels fj, Qj only and an extra coefficient c = J^^^i ({"pp +1"'" |2p+i)2 » 
which is the same for all relays. Therefore, this power allocation can be done distributively with the 
extra knowledge of one single coefficient c, which can be broadcasted by the receiver. In our case, 
every relay has a separate power constraint. This is a more practical assumption in sensor networks 
since every sensor or wireless device has its own battery power limit. The power control solutions 
of the two cases are totally different. 

If relay selection is used and only one relay is allowed to cooperate, it can be proved easily that 
we should choose the relay with the highest 

7 1/r p\ j\jj9j\ 

We call h the relay selection function since a relay with a larger hj results in a higher receive SNR. 
While all relays are allowed to cooperate, the concepts of the best relay and relay selection function 
are not clear. Since the power control problem is a coupled one, it is hard to measure how much 
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contribution a relay has. As discussed before, in network beamforming, a relay with a larger (pj does 
not necessarily use a larger power or has more contribution. But we can conclude that if 0^ > (pi, 
the fraction of power used at relay k, ak, is no less than the fraction of power used at relay /, ai. It 
is worth to mention that in network beamforming, relays with larger enough 0's use their maximal 
powers no matter what their maximal powers are. Actually, it is not hard to see that if at one time 
channels of all relays are good, every relay should use its maximum power. 

3.3 Simulation Results 

In this section, we show simulated performance of network beamforming and compare it with per- 
formance of other existing schemes. Figures |2(a)| and |2(b)| show performance of networks with 



Rayleigh fading channels and the same power constraint on the transmitter and relays. In other 
words, fi,gi are CJ\f{0, 1) and Pq = Pi = ■ ■ ■ = Pr = P. The horizontal axis of the figures 
indicates P. In Fig. |2(a)[ simulated block error rates of network beamforming with optimal power 



control are compared to those of best-relay selection, Larsson's scheme in [24] with total relay power 
P, distributed space-time coding in [10], and amplify-and-forward without power control (every re- 
lay uses its maximal power) in a 2-relay network. The information symbol s is modulated as BPSK. 
We can see that network beamforming with optimal power control outperforms all other schemes. 
It is about 0.5dB and 2dB better than Larsson's scheme and best-relay selection, respectively. With 
perfect channel knowledge at relays, it is 7dB better than Alamouti distributed space-time coding, 
which needs no channel information at relays. Amplify-and-forward with no power control only 
achieves diversity 1, distributed space-time coding achieves a diversity slightly less than two, while 
best-relay selection, network beamforming, and Larsson's scheme achieve diversity 2. Fig. |2(b)| 
shows simulated performance of a 3-relay network under different schemes. Similar diversity re- 
sults are obtained. But for the 3-relay case, network beamforming is about 1.5dB and 3.5dB better 
than Larsson's scheme and best-relay selection, respectively. 

In Fig. |3(a)[ we show performance of a 2-relay network in which Pq = Pi = P and P2 = P/2. 
That is, the transmitter and the first relay have the same power constraint while the second relay has 
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only half the power of the first relay. The channels are assumed to be Rayleigh fading channels. In 
Fig. |3(b)[ we show performance of a 2-relay network whose channels have both fading and path-loss 
effects. We assume that the distance between the first relay and the transmitter/receiver is 1, while 
the distance between the second relay and the transmitter/receiver is 2. The path-loss exponent [38] 
is assume to be 2. We also assume that the transmitter and relays have the same power constraint, 
i.e., Pq = Pi = P2 = P. In both cases, distributed space-time coding does not apply, and Larsson's 
scheme applies for the second case only. So, we compare network beamforming with best-relay 
selection and amplify-and-forward with no power control only. Performance of Larsson's scheme 
is shown in Fig. |3(b)| as well. Both figures show the superiority of network beamforming to other 
schemes. 

4 Networks with a Direct Link 

The previous section is on power control of relay networks with no DL between the transmitter and 
receiver. In this section, we discuss networks with a DL. As in [8], there are several scenarios, which 
we discuss separately. 



4.1 Direct Link During the First Step Only 

In this subsection, we consider relay networks with a DL during the first step only. This happens 
when the receiver knows that the transmitter is in vicinity and listens during the first step, while the 
transmitter is not aware of the DL or is unwilling to do the optimization because of its power and 
delay constraints. It can also happen when the transmitter is in the listening or sleeping mode during 
the second step. 

In this case, (3o = 0. From (JH) and Q, the system equations can be written as 



Xi 
X2 



aoVPofo 



s + 



W2 






Wi 
y/l+al\n\^Po 



g-jarg/. 
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Using maximum ratio combining, the ML decoding is 

2 



arg max 



xi -aoVPofos 




o/ -^ 9ll9r-)\ ^ l_fl /TT 



^^T^ 



:S 



The optimization problem is thus the maximization of the total receive SNR of both transmission 
branches, which equals 

2 






\ha^\VPi_ 

■|2Pn 

■2^1x121 2r:)\ V-' -\i'i^\ u 



First, both terms in the SNR formula increase as a^ increases. Thus, a^ = 1, i.e., the transmitter 



should use its maximum power. The SNR optimization problem becomes the one in Section [37U in 
which there is no DL. Therefore, the power control of networks with a DL during the first step only 
is exactly the same as that of networks without a DL. This result is intuitive. Since with a DL during 
the first step only, operations at both the transmitter and relays keep the same as networks without 
the DL. The only difference is that the receiver obtains some extra information from the transmitter 
during the first step, and it can use the information to improve the performance without any extra 
cost. For the single-relay case, it can be proved easily that to maximize the receive SNR, the relay 
should use its maximal power as well, that is, al = 1. 

4.2 Direct Link During the Second Step Only 

In this subsection, we consider relay networks with a DL during the second step only. This happens 
when the transmitter knows that the receiver is at vicinity and determines to do more optimization 
to allocate its power between the two transmission steps. However, the receiver is unaware of the 
DL and is not listening during the first step. It can also happen when the receiver is in transmitting 
or sleeping mode during the first step. 

In this case, a;i = and X2 is given in (|3]). The receive SNR can be calculated to be 

2 



P I R \ f \ -L n, \^^ 0'i\fi9i\'/P'i 



1 + Er=i T 
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First, we show that a^ + Pq should take its maximal value 1, i.e., the transmitter should use all its 



power. Assume that dg + /9q < 1 is the optimal solution. Define /3o = a/1 — a^. We have /?o > Pq. 
Therefore, SNR{ao, /3o) < SNR{ao, (Sq). This contradicts the assumption that (do, Po) is optimal. 
Define 

\9iWPi I ao\fi9iWPi . k 



tti = — •" • " , bi = — , Cj = ^^, A = diag-^a^ and % = a,- a. 

The receiver SNR can be calculated to be 



(yr^^i/oi + (b,x))' (yr^^i/oi + (c,y))' 



V'(ao,x) = Po ^ „ -i „, = ^c 



i + ||ixp i + l|yP 

For any fixed ao, we can optimize ai, . . . , a/? following the analysis in Section [3TT1 The following 
theorem can be proved. 

Theorem 2. Define (pj = ^far j = 1, . . . ,R and 4>r+i = 0. For any fixed a^ G (0, 1), order (pj as 

For i = 0, . . . , R, let Xi = . """^ "^'" , — an J define x*^*' ?5 defined as 



.(j) 
4 



1 j=fl,...,fi 



r/?e receive SNR is maximized at x^*"\ where Iq is the smallest i such that A.j < 0r^ . 

Proof. The proof of this theorem follows the one of Theorem [T] and the lemmas it uses. D 

As discussed in Section [BTTI for networks with no DL, there is no need to consider the solution 
of Subproblem 0. Here it is different. Define fi = 0!^^^ || c|| . If we denote the solution of Subproblem 

0, max|y|g[o,fi],Oi{dyda i+r ya ' ^^ y^^'*' because of the existence of the DL during the second step, 
it is possible that ||y^°''|| < ^i. 



Now we discuss the optimization of ao- We first consider the case of ao G (0, 1). For any given 

r 1 -^ 

= ai" ■ ■ aR J 



the ao that maximizes the receive SNR satisfies tt^ = 0. Thus, the optimal ao 



dao 
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can be found numerically by solving -^J^ = 0. It can be proved easily that ^ > when ao — * 0^ 



■^ < when ao ^ 1". Thus, the maximum of ^ is reached inside (0, 1). 
When the power at the transmitter is high (Pq ^ l)^ the receive SNR can be approximated by 

2 



ip{ao,x) ^ d{ao) = . , , , 2 

1 + d2/aQ 



where di = -h=^■^-^^ai\gi\^/Pi, and d2 = ■p-J2i=i'^'i\9i/fi\^Pi- It can be calculated straightfor- 
wardly that for tto G (0, 1), 



dd 



and 

dd 



■\fo\ao - 2&|/o|«o + ^l/ol + «^\/ 1 - «o 



dao 



^ Ifol'al - 4d2|/or«o + 2rf2|/orao + c?^(4|/or - c?I)aS + dM - Ifol') = 0. 



This is a quartic equation of «q, whose solutions can be calculated analytically. Note that -^ > 
when ao -^ 0+ and ^ < when ao -^ 1^. Thus the maximum of d is reached inside (0, 1). An 
approximate solution of ao can thus be obtained analytically at high transmit powers. 

Now we consider the cases of ao = and ao = 1- If ^o = 0, the system degrades to a point- 
to-point one since only the DL works. Thus, the receive SNR is |/opPo. For oq = 1, we can 
obtain the optimal x using Theorem [21 Thus, we obtain three sets of a and x for the three cases: 
tto £ (0, 1), ao = 0, and «o = !> respectively. The optimal solution of the system is the set of ao 
and X corresponding to the largest receive SNR. The power control problem in networks with a DL 
during the second step only can thus be solved using the following recursive algorithm. 

Algorithm 1. 

1. Initialization: Set yif^^^^°^'^ = 1/}, the R-dimensional vector of all ones, g]\j]^Y'^'"'^°'^'^i = q^ 
and count = 0. Set the maximal number of iterations iter and the threshold thre. 

2. Optimize ao with x = ■y^(P^'^'"'^°^'^>_ Denote the solution as a^ . We can either do this numeri- 
cally or calculate the high SNR approximation. 
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.(1) 



3. With ao = ttg > fi^d. the x that maximizes the receive SNR using Theorem^ Denote it as xi. 



(1) 



Calculate SNRi = ^(ag , x 



4. Set count = count+1. If count < iter and 



SNRi - SNR'^^ 



(previous) 



> thre, setx}( 



(previous) 



xi, SNR\ 



(previous) 



SNRi, and go to step 2. 



5. Find the solution of^ with ag = 1 using Theorem^ Denote this solution as X2. 

6. The optimal solution is: (og, x*) = argmax < ^^(ag , Xi), '?/;(l, X2), '?/^(0, 0/j) >. 



Similarly, the distributive strategies proposed in Section |3^ can be applied here. 

4.3 Direct Link During Both Steps 

In this subsection, we consider relay networks with a DL during both the first and the second steps. 
This happens when both the transmitter and the receiver know that they are not too far away from 
each other and decide to communicate during both steps with the help of relays during the second 
step. 

From (HI) and ([3]), the system equation can be written as 



Xi 
X2 



Wi 






ttovA/o 

Similar to the networks discussed in Section 14.11 the maximum ratio combining results in the 
following ML decoding: 



arg max 



xi - aoV-Po/os 



+ 



X2 






l/»g»lv^ 



Ifil^Po 



R o.}\g^?Pi 



1 I Y^-K "jlffil -f^^i 



The total receive SNR of both transmission branches can be calculated to be 

Pof/5o|/o|+«oEii^S^)' 

1 + Ei=l l+a^l/iPPo 
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The same as the networks in Section l4~2l a^ + /3q should take its maximal value, which is 1 . That is, 



/?o = a/1 — al- Similar to the SNR optimization in Section HTH for any given ao G (0, 1), the SNR 
maximization is the same as the maximization of ip, which is solved by Theorem [2l But due to the 
difference in the receive SNR formula, the optimal ao given ai, . . . , a^j is different. It is the solution 
of 2ao-Po|/oP + g^ = 0. When the DL exists during both steps, the case of ao = 0, whose receive 
SNR is I /o P-Po will never outperform the case of ao = 1 , whose receive SNR is \fo\'^Po + ip{l,x.){or 
some X. Thus, the case ao = needs not to be considered. The power control problem in networks 
with a DL during both steps can thus be solved using the following recursive algorithm. 

Algorithm 2. 

1. Initialization: Set x^^'"'^^*""'*^ = 1^, ^j^^j^yprevious) _ q^ ^^^ count = 0. Set the maximal 
number of iterations iter and the threshold thre. 

2. Optimize a^ with x = ^(p^'^'"'-°^'^\ Denote the solution as a^ . We can do this numerically. 

3. With ao = ag , find the x that maximizes ip using Theorem |2] Denote it as xi. Calculate 

SNR,= (^ai^^yWPo + ^iai^\^i). 



^{previous) 



(previous) 



4. Set count = count+1. If count < iter and SNRi — SNR{^ > thre, setx^^ 
xi, ^iVi?;^"^"°"'') = SNR, and go to step 2. 

5. Find the solution ofx with ao = 1 using Theorem^ Denote this solution as X2. 

6. The optimal solution is: (ap , x* ) = arg max < ( ag ) | /o P-Po + "V^ (^o , xi ) , | /o pPo + "^ ( 1 > ^2 ^ 

Again, the distributive strategies proposed in Section [X2l can be applied here. 

4.4 Performance Comparison 

In this subsection, we compare single-relay networks in which the power constraints at the transmit- 
ter and the relay are same, i.e., Pq = Pi = P. The channels are assumed to have both the fading 
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and path-loss effect. There are four cases: no DL, a DL during the first step only, a DL during the 
second step only, and a DL during both steps. 

In Fig. [51 we compare networks in which the distance of every link is the same, i.e., the three 



nodes are vertexes of an equilateral triangle with unit-length edges as shown in Fig. |4(a)[ We can see 
that the network with no DL has diversity 1 while networks with a DL and power control achieve 
diversity 2. The network with a DL during the first step performs less than 0.5dB better than the 
network with a DL during the second step only, while the network with a DL during both steps 
performs the best (about IdB better than the network with a DL during the first step only). To 
illuminate the effect of power control, we show performance of networks whose transmit power at 
the relay and transmitter are fixed. For the network with a DL during the first step only, there is no 
power control problem since it is optimal for both the transmitter and the relay to use their maximal 
powers. For the other two cases, we let the transmitter uses half of its power, P/2, to each of the two 
steps and the relay always uses its maximum power P. We can see that, if the DL only exists during 
the second step, without power control, the achievable diversity is 1. At block error rates of 10"^ 
and 10^^, it performs 3 and 6dB worse, respectively. For networks with a DL during both steps, 
power control results in a 1.5dB improvement. 



In Fig. |6(a)| and |6(b)[ we show performance of line networks with path-loss exponents 2 and 3 



respectively. As shown in Fig. |4(b)[ the three nodes are on a line and the relay is in the middle of 
the transmitter and receiver. The distance between the transmitter and receiver is assumed to be 2. 
The same phenomenon as in the equilateral triangle networks can be observed. The network with 
a DL during both steps performs the best (about IdB better than the network with a DL during the 
first step only). The network with a DL at first step only performs slightly better than the one with a 
DL during the second step only. But the difference is smaller than that in Fig. [5l The performance 
difference between line networks with and without DLs is smaller than those in equilateral triangle 
networks, and it gets even smaller for larger path-loss exponents. This is because as the distance 
between the transmitter and receiver or the path-loss exponent is larger, the quality of the DL is 
lower. Therefore, the improvement due to this link is smaller. For both cases, power control results 
in a 1.5dB improvement when the DL link exists for both steps and a higher diversity when the DL 
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exists for the second step only. 



Then we work on the random network in Fig. |4(c)[ in which the relay locates randomly and 
uniformly within a circle in the middle of the transmitter and the receiver. The distance between the 
transmitter and the receiver is assumed to be 2. The radius of the circle is denoted as r. We assume 
thatr < 1. This is a reasonable model for ad hoc wireless networks since if communications between 
two nodes is allowed to be helped by one other relay, one should choose a relay that is around the 
middle of the two nodes. In other words, the distance between the relay and the transmitter or 
receiver should be shorter than that between the transmitter and receiver. 



We work out the geometry first. As in Fig. |4(c)[ we denote the positions of the transmitter. 



the receiver, the relay, and the middle point of the transmitter and the receiver as A, C, D, and B, 
respectively. Denote the angle of AB and BD as 9 and the length of BD as p. The lengths of AD 



and CD are thus a/1 + p^ — 2pcos6 and yl + p^ + 2pcos6. Since D is uniformly distributed 
within the circle, 9 is uniform in [0, tt) and the pdf and cdf of p can be calculated to be 



2 



pip) = — and P{p < x) = — 



respectively. Define Y = ry/X. If X is uniform on (0, 1), it can be proved that 

P(y <x) = FirVx <x) = F (x <^] =^. 

Thus, Y has the same distribution as p. Therefore, we generate Y to represent p. 

Fig. |7] shows performance of random networks with path-loss exponent 2 and r = 1/2. We can 
see that the same phenomenon as in line networks can be observed. With a DL at both steps, the 
random network performs about IdB worse than the line network. 

5 Conclusions and Future Work 

In this paper, we propose the novel idea of beamforming in wireless relay networks to achieve 
both diversity and array gain. The scheme is based on a two-step amplify-and-forward protocol. 
We assume that each relay knows its own channels perfectly. Unlike previous works in network 
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diversity, the scheme developed here uses not only the channels' phase information but also their 
magnitude. Match filters are applied at the transmitter and relays during the second step to cancel the 
channel phase effect and thus form a coherent beam at the receiver, in the mean while, optimal power 
control is performed based on the channel magnitude to decide the power used at the transmitter and 
relays. The power control problem for networks with any numbers of relays and no direct link is 
solved analytically. The solution can be obtained with a complexity that is linear in the number of 
relays. The power used at a relay depends on not only its own channels nonlinearly but also all other 
channels in the network. In general, it is not even a differentiable function of channel coefficients. 
Simulation with Rayleigh fading and path-loss channels show that network beamforming achieves 
the maximum diversity while amplify-and-forward without power control achieves diversity 1 only. 
Network beamforming also outperforms other cooperative strategies. For example, it is about 4dB 
better than best-relay selection. 

Relay networks with a direct link between the transmitter and receiver are also considered in this 
paper. For networks with a direct link during the first step only, the power control at relays and the 
transmitter is exactly the same as that of networks with no direct link. For networks with a direct 
link during the second step only and networks with a direct link during both steps, the solutions are 
different. Recursive numerical algorithms for the power control at both the transmitter and relays are 
given. Simulated performance of single-relay networks with different topologies shows that optimal 
power control results in about 1 .5dB improvement in networks with a direct link at both steps and a 
higher diversity in networks with a direct link at the second step only. 

We have just scratched the surface of a brand-new area. There are a lot of ways to extend and 
generalize this work. First, it is assumed in this work that relays and sometimes the transmitter 
know their channels perfectly, which is not practical in many networks. Network beamforming 
with limited and delayed feedback from the receiver is an important issue. In multiple- antenna 
systems, beamforming with limited and delayed channel information feedback has been widely 
probed. However, beamforming in networks differs from beamforming in multiple-antenna systems 
in a couple of ways. In networks, it is difficult for relays to cooperate while in a multiple-antenna 
system, different antennas of the transmitter can cooperate fully. There are two transmission steps in 
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relay networks while only one in multiple-antenna systems, which leads to different error rate and 
capacity calculation and thus different designs. Second, the relay network probed in this paper has 
only one pair of transmitter and receiver. When there are multiple transmitter-and-receiver pairs, an 
interesting problem is how relays should allocate their powers to aid different communication tasks. 
Finally, the two-hop protocol can be generalized as well. For a given network topology, one relevant 
question is how many hops should be taken to optimize the criterion at consideration, for example, 
error rate or capacity. 
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Figure 1 : Wireless relay network. 
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(a) 2-relay network (b) 3-relay network 

Figure 2: Networks with fading channels and same power constraint for all nodes. 
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(a) Network with different relay powers (b) Network with path-loss plus fading channels 

Figure 3: 2-relay networks with different relay powers and pass-loss plus fading channels. 
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Figure 4: Network topology. 
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Figure 5: Equilateral triangle network. 
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(a) Path-loss exponent 2 (b) Path-loss exponent 3 

Figure 6: Single-relay line network. 
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Figure 7: Single-relay random network with pass-loss exponent 2 and r = 1/2. 
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